Transport signatures of pseudo-magnetic Landau levels in strained graphene ribbons 
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In inhomogeneously strained graphene, low-energy electrons experience a valley-antisymmetric 
pseudo-magnetic field which leads to the formation of localized states at the edge between the 
valence and conduction bands, understood in terms of peculiar n = pseudo-magnetic Landau 
levels. Here we show that such states can manifest themselves as an isolated quadruplet of low- 
energy conductance resonances in a suspended stretched graphene ribbon, where clamping by the 
metallic contacts results in a strong inhomogeneity of strain near the ribbon ends. 
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Graphene [T is a one- atom- thick crystalline mem- 
brane [2 capable of withstanding reversible deformations 
of up to 10% [3 , which is important because of the pecu- 
liar way strain affects the electronic properties of this ma- 
terial. It is a common feature of all materials with several 
degenerate valleys in the band structure [4] that the effect 
of their lattice deformations on electrons is equivalent to 
that of an effective gauge field [3l[5l[6]. In graphene, elec- 
trons near the Fermi level occupy states in the vicinity of 
a Dirac point at the edge between the valence and con- 
duction bands, in one of two inequivalent valleys centered 
at the corners of the hexagonal Brillouin zone, K and K' . 
Consequently, inhomogeneous strain in graphene influ- 
ences electron motion in a manner similar to an effective 
pseudo- magnetic field, which has the opposite sign in the 
two valleys [7]. Recent scanning-tunneling experiments 
on graphene nanobubbles [8 have shown that even for 
relatively weak deformations such pseudo-magnetic fields 
can reach values equivalent to tens, and even hundreds of 
Tesla, resulting in the formation of a discrete spectrum 
of 'Landau levels' (LL), including the peculiar n = LL 
state pinned to the edge between the valence and con- 
duction bands. 

Here we show that such pseudo-magnetic LL states 
form in the contact regions of stretched graphene 
nanoribbons (GNR) and then give rise to a characteristic 
signature in the electronic transport: a quadruplet of low- 
energy conductance resonance, slightly split by the valley 
mixing and the tunnel coupling via evanescent modes in 
the middle part of the GNR. 

The considered device is comprised of a GNR that is 
clamped at the ends, and contacted by heavily doped 
unstrained graphitic leads. We choose the ribbon to have 
armchair side edges along the transport direction and 
set contacts with bulk electrodes along the y direction. 
Such a ribbon can be described using the tight-binding 
Hamiltonian ^ 

'^ = E^^4^^ + E^^^-4c,. (1) 

^ {ij) 

where q is a fermionic annihilation operator acting on a 



site i and {ij) denote pairs of nearest neighbors. Com- 
pared to pristine monolayer flakes, the on-site potential 
Vi is modulated by strain, which we take into account by 
Vi = ^^divu(ri) where u = [ux^Uy) is the displace- 
ment field of the membrane and Cc is the on-site energy 
of electrons in a lattice with a given carbon-carbon bond 
length r = 1.42 A. The hopping matrix elements 

1^J = 7oe^°^^^^/""'\ kj - r(l + n,, • wn,j), (2) 

depend on the distance lij between lattice sites, modi- 
fied by the strain |3|. Here w is the 2x2 strain tensor 
'^Q!/3 = \{dcxUj3 + dpUoc) with = X or the param- 
eter 7^0 = ^mnr ^ ~3 relates the change of the nearest 
neighbor coupling to the change of the bond length [9] 
(70 - -3eV), and n,, = (1, 0), (-^, and (-^ -^) 
are the unit vectors along carbon-carbon bonds in the 
honeycomb lattice. For a homogeneously strained arm- 
chair GNR, Wxx = u)^ u)yy = —crw^ and Wxy = 0, where 
(J = 0.165 is the Poisson ratio for graphite JOj and w 
parameterizes tensile strain. The strain-induced asym- 
metry in the hoppings between neighboring carbon sites 
is equivalent to the effect of a valley-dependent vector 
potential 

where ^ = ±1 for valleys K and K' . Clamping the ends 
of a suspended strained ribbon makes the strain distri- 
bution near the device contacts inhomogeneous [H], and 
this results in the appearance of a pseudo-magnetic field 
3 = rotv4 in the vicinity of the contacts [71 116] . This is 
illustrated in Fig. [TJa), where we applied linear elasticity 
theory [151 HI] to three ribbons of width W 40nm and 
aspect ratio L/W = 2, 3, and 4, with w = 0.05 in their 
middle parts. The pseudo-magnetic field is the largest 
positive (blue) or negative (red) near the contacts at the 
right and left ends, and is small in the middle of the 
ribbon. 

The phase-coherent transport properties of two- 
terminal devices are encoded in the scattering matrix [17] 
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FIG. 1. We consider the transport through strained suspended graphene nanoribbons (GNR) which are clamped at highly- 
doped contacts. Panel (a) shows the distribution of pseudo-magnetic fields B{T) for electrons in the K valley for GNR with 
W ^ 40nm and aspect ratios L/W = 2, 3, and 4. The inhomogeneous tensile strain in the middle of the nanoribbon is 
w = 0.05. Panel (b) shows the spatial structure of electron wave amplitudes corresponding to several resonances identified in 
panel (d), which displays the zero-temperature conductance of the ribbons as a function of the Fermi energy. For comparison, 
panel (c) shows the conductance for the ribbon with L/W = 3 and no strain (w = 0, blue) or artificially imposed homogeneous 
strain {w = 0.05, red). 



which we evaluate using the recursive Green's function 
technique ^3 [19], applied to the tight-binding model 
of the GNR sketched in Fig. [2] Within the Landauer- 



Biittiker formalism [20] we calculate the device conduc- 
tance as a function of the Fermi level in the middle 
part of the GNR, for a given height {V = -200meV) 
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FIG. 2. Sketch of the tight-binding model ^ of the GNR 
junction with armchair boundaries. The system is comprised 
of two ideal heavily doped leads {V = — 200meV) and a cen- 
tral suspended region, in which strain modulates the hopping 
matrix elements jij and the on-site energy Vi. 

of the gate-controlled potential energy step between the 
doped graphene leads and the suspended part. The 
resulting device is an n-p-n {Ep < OmeV) or n-n'-n 
{Ef > OmeV) graphene junction, where most of the con- 
ductance features are determined by scattering from the 
strain-modified n-p or n-n^ interfaces. 

Figure [TJc) shows the Fermi energy dependence of 
the zero-temperature conductance of the ribbon with 
= 3, for no strain (blue) and for artificially imposed 
homogeneous strain with w = 0.05 (red). The unstrained 
armchair GNR is semiconducting; the conductance there- 
fore exhibits a gap around the DP {Ef = OmeV). The 
conductance oscillations away from the DP are due to the 
Fabry-Perot-like standing wave resonances in the electron 
transmission across the potential barrier geometry [21 . 
For w = 0.05 homogeneous strain the conductance is 
completely suppressed for \Ef\ < lOOmeV. This is be- 
cause the constant vector potential induced by the ho- 
mogeneous strain shifts the Dirac cones away from the 
K and corners of the Brillouin zone, perpendicular to 
the transport direction [3 . The misalignment between 
the Fermi surfaces in the unstrained leads and in the 
strained suspended region results in a suppression of the 
conductance in the ballistic regime [22 . The threshold 
for such an insulating behavior for parameters used in 
Fig. [TJc) is w = 0.024, and can be lowered by reducing 
the height of the potential step V between the central 
part of the ribbon and the contacts. 

This behavior is significantly changed when the in- 
homogeneity of the strain is taken into account. Fig- 
ure [ijd) shows the Fermi-energy dependence of the 
conductance for the inhomogeneously strained ribbons 
shown in Fig. [ija). In contrast to Fig. [ijc), here, we 



FIG. 3. Sublattice-resolved electron amplitude for one of the 
resonances m Fig.[l];d) {L/W = 2,E = -5.84meV), obtained 
by placing the probing perturbation on the A (left panel) or 
B (right panel) sites. 

find several groups of additional sharp resonance con- 
ductance peaks in the energy range \Ef\ < lOOmeV. To 
reveal the nature of each group of these peaks, we an- 
alyze the spatial distribution of the corresponding elec- 
tronic states. Within the Landauer-Biittiker formalism, 
this can be obtained from the response to local perturba- 
tions of the scattering amplitudes at the energies close to 
the resonance conditions [S]- The reconstructed spatial 
structure of the resonance states is shown in Fig. [TJb) . 
As illustrated in the top two rows, the states away from 
the DP correspond to Fabry-Perot-like standing waves 
formed in the homogeneously strained central part of the 
structure, to which they are confined due to electron re- 
flection from the interfaces. These resonances appear be- 
cause the transverse momentum is no longer conserved 
when the interface region is inhomogeneous, which al- 
lows carriers to overcome the misalignment between the 
Fermi surfaces described above. However, two groups of 
resonances, in the energy range — 25meV< Ef < OmeV 
just below the DP, display a very different behavior. As 
shown in the insets in Fig. [ijd), these resonances occur 
in almost degenerate pairs, of which we find two groups. 
The spatial structure of these states, shown in the bot- 
tom four rows of Fig. [TJb), clearly resembles the pseudo- 
magnetic field distribution. 

We now demonstrate that this quadruplet of reso- 
nances can be attributed to the n = pseudo-magnetic 
Landau level. Our main piece of evidence is a unique 
feature of this LL in graphene, namely that the elec- 
tron amplitude resides either on the A or B sublattice. 
The selected sublattice depends on the sign of B but is 
independent of the valley |5j. In contrast, higher LLs 
and Fabry-Perot-like resonances occupy both sublattices 
equally. By selectively placing the probing perturbation 
on the A OT B sites, we find, as illustrated for one ex- 
ample in Fig. |3j that the low-energy resonances have the 
following property: their amplitude is high on the A-sites 
near the left end of the ribbon (where S < 0), and the 
amplitude is high on the B-site near the right end (where 
B > 0). This fully agrees with the unique feature of the 
n = LL described above. 

The fact that we find four such low-energy resonances, 
as well as the dependence of their splitting on the length 
of the ribbon, further supports this interpretation of the 
origin of these states. The y —y reflection symmetry 
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maps valleys K and K' onto each other. The resulting 
symmetric and anti-symmetric superpositions of the two 
valley manifestations of the n = LL lead to a splitting of 
the quadruplet into two groups: a low-energy branch at 
Ep ~ — 24meV, which is valley-symmetric and displays 
a maximum on the symmetry axis, and a high-energy 
branch dX Ep ^ — 7meV, which displays a nodal line on 
this axis. Each group splits further into two narrowly 
spaced lines because of the tunnel coupling of the states 
and near the left and right ends of the ribbon, 
which is provided by the evanescent tails of the electronic 
wave functions in the middle part, where B is small. Note 
that the splitting in each of these pairs is smaller in a 
longer ribbon, since the overlap of the evanescent tails of 
^ L and ^ R decreases with the separation between the 
GNR ends. 

In conclusion, we describe a unique transport signa- 
ture of the pseudo-magnetic field in a strained suspended 
graphene nanoribbon, namely the resonant transmission 
via the sublattice-polarized n = pseudo-magnetic Lan- 
dau level. These states form near the inhomogeneously 
strained contact regions and give rise to a characteris- 
tic quadruplet of conductance resonances near the Dirac 
point. The above-proposed analysis is directly appli- 
cable to graphene ribbons where high-quality armchair 
edges are obtained by the oriented growth on patterned 
SiC substrates [24], etching of graphene samples with 
catalytic nanoparticles [25^, or chemical derivation [26] . 
Even though imperfections in the system will lead to 
the appearance of additional conductance resonances, the 
n = Landau level state is protected against the infiu- 
ence of disorder by its unique position at the energy of the 
Dirac point and its energetic separation from the Fabry- 
Perot-like resonances. 
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